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Abstract 

The ^5o pairing gap in isospin-symmetric nuclear matter and finite nuclei is investigated using 
the chiral nucleon-nucleon potential at the N'^LO order in the two-body sector, and the N^LO 
order in the three-body sector. To include realistic nuclear forces in RHB (Relativistic Hartree 
Bolgoliubov) calculations we rely on a separable form of the pairing interaction adjusted to the 
bare nuclear force. The separable pairing force is applied to the analysis of pairing properties for 
several isotopic and isotonic chains of spherical nuclei. 
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I. INTRODUCTION 



Nuclear energy density functionals (EDF) provide a microscopic, globally accurate, and 
yet economic description of ground-state properties and collective excitations over the whole 
nuclide chart. Even though it originates in the effective interaction between nucleons, a 
generic density functional is not necessarily related to any given nucleon-nucleon (NN) po- 
tential and, in fact, some of the most successful modern functionals are entirely empirical 
|l| . Of course it is very much desirable to have a fully microscopic foundation for a universal 
density functional, and this is certainly one of the major challenges in low-energy nuclear 
structure physics |2|. 

The EDF approach to nuclear structure is analogous to Kohn-Sham density functional 
theory (DFT) and, because it includes correlations, goes beyond the Hartree-Fock approx- 
imation. Kohn-Sham DFT has the advantage of being a local scheme, but its usefulness 
crucially depends on our ability to construct accurate approximations for the most impor- 



In a series of recent articles 



tant part of the functional, that is, the universal exchange-correlation functional js] 



6| concepts of effective field theory and density functional 
theory have been used to derive a microscopic relativistic EDF-based model of nuclear many- 
body dynamics constrained by in-medium QCD sum rules and chiral symmetry. The density 
dependence of the effective nucleon-nucleon couplings in this model (called FKVW in the fol- 
lowing) is determined from the long- and intermediate-range interactions generated by one- 
and two-pion exchange processes. They are computed using in-medium chiral perturbation 
theory, explicitly including A(1232) degrees of freedom [3]. Divergent contributions to the 
nuclear matter energy density, calculated at three-loop level, are absorbed by few contact 
terms. These constants are understood to encode unresolved short- distance dynamics. 

The relativistic FKVW model has successfully been employed in studies of ground-state 
properties of spherical and deformed nuclei using the Relativistic Hartree-Bogoliubov frame- 
work (RHB 8|). In the description of open-shell nuclei, in particular, a hybrid model has 
been used with the FKVW Kohn-Sham potential in the particle-hole (ph) channel and, like 
in most applications of RHB-based models, the pairing part of the empirical Gogny force 9| 
in the particle-particle (pp) channel. 

Even though this approach has been very successful, it is not theoretically consistent 
because of the choice of the empirical effective interaction in the pp channel. Quite recently, 



as part of a larger program to develop a framework of fully microscopic nuclear energy density 
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functionals, much effort has been devoted to designing non-empirical pairing functionals 

Q. 

The aim of this work is to formulate a consistent microscopic framework for open-shell 
nuclei, in which both the ph and the pp channels of the effective inter-nucleon interaction are 
determined by chiral pion-nucleon dynamics. Thus we consider a separable pp interaction 
based on a microscopic pairing interaction constrained by chiral dynamics (see Ref. 15 1 
for previous calculations involving the N^LO chiral potential), combine it with the FKVW 
functional in the ph channel and, employing the corresponding RHB model, present a study 
of pairing gaps in isotopic and isotonic chains of spherical open-shell nuclei. 

We will use the realistic NN potential developed by the Idaho group at next-to-next-to- 
next-to-leading order (N^LO) in the chiral expansion 16| (see also Ref. 17|), and a two-body 
density- dependent potential derived from the relevant diagrams at the N^LO order in the 

22 | for similar approaches and pertinent details). 



18| (see also Refs. 



Q,[l9[j22 



three body sector 

The paper is organized as follows. In Sec. Hi] we discuss results for the pairing gap of 
nuclear matter in the BCS approximation. Sec. Illll briefiy reviews the method introduced by 
Y. Tian et al. 23, 2^ to apply realistic pairing interactions to calculations of finite nuclei. 
In Sec. |IV]we analyze pairing gaps in spherical nuclei for several isotopic and isotonic chains. 
Sec. |V] summarizes the principal results. 



II. PAIRING GAP IN A HOMOGENEOUS INFINITE SYSTEM 



The momentum and density-dependent pairing field A(/c, kp) in infinite matter is deter- 
mined by the solution of the BCS gap equation 

1 r p'^V{p,k)A{p,kF) 



dp 



(1) 



^[£{p, kp) - S{kF, kp)f + A{p, kpf 

where V{p^ k) represents the off-shell pairing potential in momentum space, £'(p, kp) is the 
quasiparticle energy, and £{kp, kp) is the Fermi energy. 

The effective force in the pairing channel is, in principle, generated by the sum of all 
particle-particle irreducible Feynman diagrams 25|]. In most application to nuclear and 
neutron matter, however, only t^ lowest-order term, that corresponds to the bare nucleon- 

Terms of higher order in the effective pairing interaction 



nucleon interaction, is retained 



represent 



Refs. 
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screening corrections to the bare force, caused by medium polarization effects (see 



and references therein). In the present analysis we only consider the bare 
interaction, while a study of polarization effects will be carried out in a forthcoming paper. 
For the pairing potential V{p, k) we employ the simple ansatz: 

V{p, k) = V2b{p, k) + ViB{v, k,m) ~ V^b^V, k) + V"2ij(fcF,P, k) , (2) 

where the three-body potential is approximated by an effective two-body density- dependent 
potential V2B derived by Holt et al. in Ref isl. These authors showed that in the singlet 
channel (^^o) the overall effect of V2B{kp,p, k) is to reduce the strong S-wave attraction 



cf. Fig. 6 of Ref. [18[). As suggested in Ref. 



here we neglect a possible isotopic 



dependence that, in any case, is expected to be rather small for the nuclei considered in the 
present analysis (see also Ref. [l^). For both terms in Eq. (|2]) we follow standard procedures 
for the regulator functions, and refer the reader to the original papers for details. 

For the single-particle spectrum that appears in the denominator of the gap equation ([T]) 
we employ the simple quadratic form 

s{p, kp) - s{kF, kp) = f~!:'\ . (3) 



2M*{kF) ■ 

This approximation should suffice because the momenta p around kp give the dominant 



contribution to the integral in Eq. ([T]). 
a very recent calculation (Fig. 9 of Ref. 



'he effective nucleon mass M*{kp) was obtained in 



29|). in which the nuclear energy density functional 



was derived to first order in the two- and three-nucleon interaction using a density matrix 



expansioi 

Fig. [1] displays the pairing gap Alkpjkp) in symmetric nuclear matter as function of 
the Fermi momentum kp. We plot results of the complete calculation that includes two 
and three-body forces (solid curve), and the pairing gap obtained with only the two-body 
NN potential at N^LO (dashed curve). Our results are shown in comparison with those 



obtained in Ref. 



32| using the Viowk potential (with single particle energies computed in 



Brueckner-Hartree-Fock theory). 

— n 

^ In Ref. 29[ the two-body interaetion comprises long-range one- and two-pion exchange contributions 
and a set of contact terms contributing up to fourth power in momenta (N'^LOW potential developed 
lowering the cut-off scale to A = 414 MeV). In addition, the authors employ the leading order chiral 
three-nucleon interaction with the corresponding parameters ce, cd and ci.3_4 adjusted in calculations of 
few-body systems. Even though the results are in good agreement with previous calculations 3^, one 



should note that higher-order corrections could have non-negligible effects 
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III. MAPPING PROCEDURE 



le pairing channel of the RHB frame- 



23l 



24j where a separable form of the 



To implement the chiral NN potential at N^LO in t 
work for finite nuclei, we adopt the approach of Refs. 
pairing interaction was introduced, with parameters adjusted to reproduce the pairing prop- 
erties of the Gogny force in nuclear matter. In nuclear matter the pairing force is separable 
in momentum space: 

{k\ V'^'' \k') = -Gp{k)p{k') . (4) 

By assuming a simple Gaussian ansatz p{k) = e"*^^^^, the two parameters G and a have been 
adjusted to reproduce the density dependence of the gap at the Fermi surface, calculated 
with a Gogny force. For the DIS parameterization 9|] of the Gogny force: G = 728 MeVfm^ 
and a = 0.644 fm. Here we apply the same procedure to the chiral NN potential at the N^LO 
order. For finite nuclei, when the pairing force Eq. (jl]) is transformed from momentum to 
coordinate space, it takes the form: 

V{n, r2, r[, r'2) =G5iR- R') P(r)P(rOi (1 - P'^) , (5) 

where R = ^ (ri + 7^2) and r = ri — r2 denote the center-of-mass and the relative coordi- 
nates, P(r) is the Fourier transform of p{k): 

and the factor 1/2 (1 — P'^) projects on the ^Sq channel. The pairing force has a finite range 
and, because of the presence of the factor 6 {R — R'), it preserves translational invariance. 
Even though 6 {R — R') implies that this force is not completely separable in coordinate 
space, the corresponding ant i- symmetrized pp matrix elements can be represented as a sum 
of a finite number of separable terms, using a method developed by Talmi and Moshinsky. 
When the nucleon wave functions are expanded in a harmonic oscillator basis, spherical 
or deformed, the sum converges relatively quickly. A relatively small number of separable 
terms reproduces with high accuracy the result of a calculation performed in the complete 



basis. We refer the reader to 23|] for more details. 

The parameters of the separable pairing force take the values: G = 892.0 MeVfm^ and 
a = 0.74 fm for the N^LO potential, and G = 1045.0 MeVfm^ and a = 0.86 fm for the 
complete potential V{p, k). We note that recently a similar approach was employed in Refs. 



m 

, where a low-rank separable representation was used to reproduce directly l^owk and 
VsAT in the ^5*0 channel (for V^n the density dependence was parametrized by a polynomial 
in the Fermi momentum). 

IV. RESULTS FOR FINITE NUCLEI 



Employing the RHB model with the FKVW functional in the ph channel and the sepa- 
rable pairing force Eq. ([5]) in the pp channel, we have calculated the self-consistent ground- 
state solutions for several sequences of isotopes (nickel, tin and lead) and isotones (A^ = 28, 

= 50 and = 82). The total binding energies and average pairing gaps are compared to 
available data in Figs. [2] and 3] The experimental masses are from Ref. 33j, and the average 



proton and neutron gaps 



34| 



A 



(7) 



are compared to empirical values determined using the 5-point formula 35| for even-even 
nuclei 

A(5)(iVo) = -l [E{No + 2) - 4E{No + 1) + QE{No) - 4E{No - 1) + E{No - 2)] . (8) 

o 

E{Nq) denotes the experimental binding energy of a nucleus with A^o neutrons {Zq for 
protons). In Eq. ([7]) the sum is over proton or neutron canonical states, A^ is the diagonal 
matrix element of the pairing field in the canonical state k, and Vk denotes the corresponding 
eigenvalue of the one-body density matrix (occupation factor)^. 

The theoretical gaps shown in Figs. |5] and [3] have been calculated using the values of 
the parameters G and a that correspond to the nuclear matter pairing gaps in Fig. [H 
The gaps calculated by including only the interaction V2B (blue diamonds) reproduce on a 
quantitative level the isotopic and isotopic trends of the empirical gaps. Including the three- 
nucleon interaction V^b induces a sizable reduction of the calculat ed g aps (green diamonds). 



14l |. The calculated gaps 



and we note that similar conclusions were also reached in Ref. 
for the isotopic chains Z = 28, Z = 50 and Z = 82 indicate that missing contributions 

n n 

like, for instance particle- vibration coupling, could play an important role [36|, |37|. Fig. [3] 
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14| a somewhat different prescription was used for the pairing gaps. The theoretical gap 



2 In Rcfs. 

^Lcs (Lowest Canonical State) corresponds to the diagonal pairing matrix element in the canonical 
single-particle state whose quasi-particle energy is the lowest, whereas experimental gaps are deduced 
from binding energies using three-point mass differences centered on odd-mass nuclei. 

6 



displays similar results for the proton pairing gaps of the isotonic chains = 28, = 50 
and = 82 (we note that here the contribution of the Coulomb interaction in the pairing 



channel is neglected). The subshell c 
and at Z = 58 in the = 82 chain 



osures that appear at A^ = 40 in the nickel chain 



39|, lead to a strong reduction of pairing correlations 



in the corresponding calculated ground-states. 

The influence of three-body forces on the total binding energy is much less pronounced, 
as shown in the lower panels of Figs. |2] and [3l where we display the absolute deviations of the 
calculated binding energies from the experimental values. On the one hand this is because 
pairing correlations contribute much less than the mean-field self-energies to the total bind- 
ing. On the other hand this is a well known characteristic of a self-consistent calculation in 
that, for a given nucleus, a reduction of pairing results in an effective enhancement of the 
mean-field contribution to the total energy, and vice versa. In general, the combination of 
the FKVW ph effective interaction and the separable pairing force Eq. (jS]) produces results 
for the total binding energies that are comparable to those obtained with the best empirical 
non-relativistic and relativistic energy density functionals. For the nickel isotopes the largest 
deviations are in the region Z ^ N, where one expects additional contributions from proton- 
neutron correlations that are not explicitly included in the FKVW functional. In the tin 
isotopes the calculated masses start deviating from data in neutron-rich nuclei beyond the 
major shell closure at A^ = 82, whereas for lead nuclei the deviations are most pronounced 
in the lightest, neutron-deficient isotopes that are characterized by soft potentials and shape 
coexistence. 



V. CONCLUSIONS 



A consistent microscopic approach to the structure of open-shell nuclei has been intro- 
duced, in which both the ph and the pp channels of the effective nuclear interaction are 
fully determined by chiral pion-nucleon dynamics. By employing an ansatz for the pairing 
force that is separable in momentum space, we have performed an efficient mapping of the 
chiral potential in the pairing channel (at the N'^LO order and the N^LO in the two-body 
and three-body sectors, respectively) to an effective pp interaction for finite nuclei. The two 
parameters of the separable pairing force are adjusted to reproduce the density dependence 
of the pairing gaps in symmetric nuclear matter. The resulting effective pairing interaction 



7 



thus enables, on the one hand, the treatment of pairing correlations in finite nuclei using 
pairing functionals constrained by chiral dynamics and, on the other hand, calculations in 
the pp channel with a finite-range interaction. The significant advantage is that the com- 
putational cost is greatly reduced when compared to nonlocal finite-range forces like, for 
instance, the empirical Gogny force. A noteworthy result of the present investigation is 
that it confirms the important role of three-body forces in determining pairing gaps in finite 
nuclei. 

This work was partly supported by INFN, MIUR and MZOS (project 1191005-1010). We 
acknowledge useful discussions with N. Kaiser, W. Weise and E. Vigezzi, and would like to 
thank J. W. Holt for providing numerical values for the three-body potential. 
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FIG. 1. Pairing gaps in symmetric nuclear matter in the ^5o channel as functions of the Fermi 
momentum. The dashed curve is obtained by including only the two-body (V2_b) chiral interaction, 
whereas the solid curve also includes the contribution of three-body (V^b) forces. The gaps are 



compared to those obtained in Ref. 



32l | using the Viowk potential. 
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V^g (N3L0) 
V,,(N3LO) + V33 




FIG. 2. Theoretical average neutron pairing gaps ([7]) of the even-even isotopes of nickel Z = 28, 
tin Z = 50, and lead Z = 82, compared to the empirical values calculated from experimental 
masses using Eq. ([8]) (upper panel). Absolute deviations of the calculated binding energies from 



the experimental values 



33( 1 (lower panel). 
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